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ABSTRACT 

We present a Comptonization model for the observed properties of the energy-dependent soft/hard time lags 
and pulsed fraction (amplitude) associated with the pulsed emission of a neutron star (NS). We account for 
the soft lags by downscattering of hard X-ray photons in the relatively cold plasma of the disk or NS surface. 
A fraction of the soft X-ray photons coming from the disk or NS surface upscatter off hot electrons in the 
accretion column. This effect leads to hard lags as a result of thermal Comptonization of the soft photons. This 
model reproduces the observed soft and hard lags due to the down- and upscattered radiation as a function of 
the electron number densities of the reflector, n T f, and the accretion column, n^ ot . In the case of the accretion- 
powered millisecond pulsars IGR J00291+5934, XTE J1751-305, and SAX J1808.4-3658 the observed time 

lags agree well with the model. Soft lags are observed only if n" n 'l°' ■ Scattering of the pulsed emission in 
the NS environment may account for the observed time lags as a non-monotonic function of energy. The time 
lag measurements can be used as a probe of the innermost parts of the NS and accretion disk. We determine 
the upper and lower limits of the density variation in this region using the observed time lags. The observed 
energy-dependent pulsed amplitude allows us to infer a variation of the Thomson optical depth of the Compton 
cloud in which the accretion column is embedded. 

Subject headings: binaries: close - pulsars: individual (IGR J00291+5934, SAX J1808.4-3658, XTE J1751- 
305) - stars: neutron - X-ray: binaries 



1, INTRODUCTION 

Accreting neutron star and millisecond pulsar (MSP) 
source spectra are successfully fitted by a two-component 
model consisting of Comptonization of higher temper- 
ature blackbody photons and thermal soft X-ray emis- 
sion. The emergent radiation is presumably produced 
by thermal and dynamical Comptonization of seed pho- 
tons coming from the NS surface and an accretion disk 
[for NS: e.g. (iTitarchuk, Mastichi adis & Kvlafil 1 19961 
119971: ITitarchuk . & S haposhni kovl 120051: iPaizis et al l2006h 
and MSP: e.g. JGilfanov et alj|1998ng ierlinski et all I 2002! 
IGierlinski & Poutane n 2005; Falan ga et al.ll2005allbh l. 

Additional information for the X-ray production processes 
and emission environment can be obtained by studying the 
pulse profile and phase shift between X-ray pulses at dif- 
ferent energy ranges. For SAX J1808. 4-3658 and later for 
XTE J 175 1-305 it was found that the low-energy pulses lag 
behind the high-energy pulses (soft phase/time lags) mono- 
tonically increasing with energy and saturating at about 10- 
20 ke V (ICui etalJll998t iFordBOOOl IGierlinski & Poutanen! 
120051) . 

This time lag effect was first interpreted as a result of pho- 
ton delay due to downscattering of hard X-ray phot ons in the 
relatively cold plasma of the disk or NS surface dCui et al.l 
119981: ITitarchuk etaiH l2002h . It was argued that the pho- 
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ton time lags were an intrinsic signature of the interaction 
of the Comptonized radiation with the NS and accretion disk 
plasma. Moreover, the absolute values of the time lags (about 
hundreds yus) are consistent with the electron scattering time 
scale fc = tj(L/c). Effective Thomson optical depth of the 
cold reflector tt is of the order of a few. Typical sizes of 
the NS photosphere and half-width o f the disk L are of the 
order of 10 6 cm. On the other hand, Pouta nen & Gierlinskil 
(120031) suggested that the lags may be produced by a com- 
bination of different angular distributions of the radiation 
components and relativistic effects. Recently, soft lags were 
also found up to ~ 100 keV in the 1.67 ms accreting MSP 
IGR J0029 1+5934, where a more com plex energy depen- 
dence was revealed (Falanga et al. 20053)- 

In this Letter, we present a model for the energy-dependent 
time lags as a result of photon downscattering and upscatter- 
ing in a Comptonized medium. We show that the resulting 
time lag signs and behavior as a function of energy depend on 
the ratio of densities of the Compton emission area, NS sur- 
face area and accretion disk regions. By time lag, we mean 
the time or phase shifts between X-ray pulses at different en- 
ergies. 

2. COMPTONIZATION MODEL OF ENERGY-DEPENDENT TIME 

LAGS 

In Figure Q] we present a typical scenario of the formation 
of the X-ray spectrum of a NS binary. We assume that the 
hard radiation is produced in an accretion column or Comp- 
ton cloud (Comptonization region) and illuminates the rela- 
tively cold material of the NS and disk. The emergent hard 
X-ray emission is a result of the upscattering of NS and disk 
soft photons off of hot electrons in the Comptonization region 
and the reflection of this upscattered radiation from the NS 
and disk surfaces. The physical parameters of the source, the 
geometrical and optical thickness of the disk, the disk plasma 
conductivity and the NS magnetic field determine the posi- 
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Fig. 1 . — This cartoon illustrates the different emission patterns responsi- 
ble for the time lags of the pulsed emission where (^ ot and 6jj cf are the di- 
mensionless temperatures of the accretion column and reflector, respectively. 
Cm ~ 0.1 is the disk illumination fraction. Soft time lags of the pulsed emis- 
sion are the result of downscattering of hard X-ray photons in the relativly 
cold plasma of the disk. The fraction of hard X-ray photons ( 1 - Cm ) directly 
seen by the observer is related to the upscattered soft photons coming from 
the NS and disk. We suggest that the emergent pulsed flux (pulsed fraction) 
is formed in the Compton cloud. 

tion of the inner edge of the disk and the illumination cover- 
ing factor Cm- The fraction of the hard photons intercepted 
by the disk and NS and ultimately reflected (Cm) can be in 
the range of at most 10-15%, depending on the details of 
the geometry. The parameter C u is a product of the albedo, 
A ~ 0.3, of the cold material (Basko, Sun vaev & Titarchukl 
1 19741) and the fraction of the Comptonized photons, f mt , emit- 
ted by the isotropically radiated central source and then inter- 
cepted by the infinite disk. The fraction f mt is less then 25% 
(lLapidus.et al.lll985l) . Thus Cm = A ■ f int l0.l. In Sec. [3] 
we use Cm = 0. 1 as a representative value for low magnetized 
(~ 10 8 G) neutron stars where the innermost edge of the disk 
is just a few kilometers from the NS surface. 

The photon energy change per electron scattering SE/6u de- 
pends on the photon energy, E, and electron plasma tempera- 
ture kT e : 

±^=40-z, (1) 
Z ou 



where z = E/m e c is the dimensionless photon energy and = 
kT e /m e c 2 is the dimensionless plasma temperature with re- 
spect to the electron rest energy m e c 2 d Sunvaev & Titarchukl 
1980). In a cold plasma (NS photosphere or accretion disk) 
for which the plasma temperature kT e is much less than a 
given energy of the photon E, the solution of the differential 
equation (Q} is 

1 1-40/z. 

u — — In (2) 

40 1 - 40 /z 

where z* is the dimensionless energy of a given photon after 
u scatterings for which the initial dimensionless energy is z, 
u = (n e <Tic)t is the dimensionless time measured in units of 
mean time between scatterings (number of scatterings), n e is 
the plasma number density, and <x T is the Thomson electron 
cross-section. In a hot plasma, e.g. the Comptonization re- 
gion, kT e , is much higher than the seed photon energy E s . In 
this case, the solution of Eq. (HJ is 



u = ln(z/z s )/40. 



(3) 



The appropriate solution of the exact kinetic equation 
dTitar chuk & Lyubarskij 1 19951) gives 



In z/Zs 



where a is the energy spectral index. One can approximate 
Eq. © by Eq. (f3]l if it is assumed that (3 + a)S«l and all. 
Note that Eq. (Q]l is obtained using a diffusion approximation. 

The resulting time lags are a linear combination of the pos- 
itive (hard) time lags formed in the Comptonization emission 
area and the negative (soft) ones formed in the accretion disk 
and NS photosphere as a result of reflection of the hard radia- 
tion (see Fig. [TJi: 

Cn 
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1 -40 ref /z 
40 iei ~ 1 -4ftrf/z, 
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Cm ln[l +(3 + a)8 hot ] 



(5) 
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where rajr is the electron number density of the reflector, «^ ot 
is the electron number density of the Comptonization emis- 
sion area (accretion column) and K f = kT r e ef / (m e c 2 ) is the 
dimensionless temperature of the reflector. We assume a typ- 
ical value of ref < 0.7 keV/511 keV. kT* ot and a are the 
best-fit parameters for the hot plasma temperature and spec- 
tral index of the Comptonization spectrum respectively, and 
Shot = kT*° l /(m e c 2 ). 

We also assume that the seed photon energy is near the low- 
est energy of the downscattered photons, i.e. z s ~ Z* (see Eqs. 
[21[4]i. The observational value of E* is about 3 keV. We set the 
value of the illumination factor Cm at about 0. 1 (see above for 
details). 

The essential model parameters for fitting the observed lags 
are the number densities n x f and n^° l . Eq. (0 is insensitive 
to 6> re f if re f < z*/4. The first term in parentheses of Eq. 
§5§ is related to the downscattering lags, and the second term 
is related to the upscattering lags. The positive (hard) lags 
cannot be observed if the density of the hot plasma is much 
higher than that of the "cold" reflector. This is because the 
mean free path between consecutive scatterings / = 1 /(n c <x T ) 
and therefore the time lag Af = ul/c are inversely proportional 
to the number density n e . Consequently, the photons spend 
less time Af in the higher density plasma than in the lower 
density plasma to undergo the same number of scatterings. 

The observer detects the negative (soft) lags only if <sc 
n h e °' and thus the relative fraction of the negative (downscat- 
tering) time lags in the apparent time lag sum is 100% (i.e. 
Cm = 1). In this case, the observed absolute values of the soft 
time lags monotonically increase with energy. Using all of the 
conditions described above Eq. [5] can be written, after some 
algebra, as 

A--^r( i - i ). <B 

a-Tn^cyz, z) 

This type of time lag correlation with energy was found for 
the first time in the accreting MSP SAX J1808.4-3658 (see 
ICui et al.lfT998h . 

3. APPLICATION OF THE MODEL TO OBSERVED TIME LAGS 

To apply the Comptonization model Eq. © to ob- 
served time lags, we used the observed pulse phase lag 
data and best-fit spectral parameters kT^ ot and a of the 
accreting MSPs IGR J00291+5934, XTE J1751-305, and 
SAX J1808.4-3658. Figure |2] shows, from top to bottom, 
IGR J0029 1+5934, XTE J1751-305, and SAX J1808.4-3658 
( Falan ga et alj|2005at iGierlihski & Poutanenll200l ICui et alj 
119981; iFordl 12000). The measurements were made either 
through computing cross-power spectra between different en- 
ergy bands or through cross-correlating the folded pulse pro- 
file in different energy bands. 
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Fig. 2. — The measured soft time lags of the pulse profile versus energy 
(crosses) with respect to the first energy channel. The best-fit curve using the 
Comptonization model (see Sec. |5J is shown with the solid line. The dot- 
dashed lines in panel (a) correspond to the upper and lower limits of the elec- 
tron number densities of the Comptonization emission area, n^ ot and of the 
reflector, n' cf in IGR J0029 1+5934. The panels (b) and (c) show XTE J1751- 
305 and SAX J1808. 4-3658, respectively. The dot-dashed lines correspond 
to the upper and lower limits of «™ f . 

Figure [2] shows the best-fit curves using the Comptoniza- 
tion model described in Sec. [2] We used kTf - 0.4 keV 
as the reflector temperature for all fits. The model for fitting 
the time lags in XTE J1751-305 and SAX J1808.4-3658 has 
only one free parameter, the number density of the "cold" re- 
flector nf (see Sec. 0. Presumably, in XTE J175 1-305 and 
SAX J1808.4-3658 the density of the Comptonization region 
is much higher than that of the "cold" reflector, i.e. nf » 
nf . This may indicate that the seed photons are Comptonized 
(upscatt ered) in the very dense plasma of the accretion c olumn 
(see e.g. lBasko & Sunvaevll976HBecker & Wolffll2006l) . The 
best-fit values are nf = 6.3 x 10 19 cm" 3 and 3.2 x 10 19 cm" 3 
for XTE J175 1 -305 and SAX J 1 808.4-3658, respectively. The 
fits of the time lag data for IGR J0029 1+5934 provide us 
with best-fit values for the "cold" and hot plasma densities. 
Both the positive (upscattering) and negative (downscatter- 
ing) time lags contribute to the apparent time lag because 
nf = 6.9 x 10 18 cm" 3 and nf = 2.1 x 10 18 cm" 3 are of 
the same order of magnitude. 

The time lag data of IGR J00291+5934, XTE J1751-305, 
and SAX J1808.4-3658 were col lected in different time inter- 
vals lasting from hours to days dCui et al.lf l998; Ford] 120001; 
Gierlinski & Poutanen 2005; Falanga et al. 2005a). However, 
a hydrodynamical (density perturbation) time, thydro, m the 
innermost part of the X-ray NS source is of the order of the 
ratio of the NS radius to the sound speed, that is, f/, y£ / ra ~ 
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Fig. 3. — The observed energy dependent pulsed fraction of (a) 
IGR 100291+5934 , (b) XTE J1751-305, and (c) SAX J1808.4-3658 along 
with the best-fit model (see Sec. 4). 

^Ns/csound^ 0.1 s. Thus, during the data collection periods, 
the densities of the surrounding plasma can vary. As a re- 
sult, the time lags also vary since they are sensitive to density 
variations (see Eq. [5). Therefore, we can infer these den- 
sity variations from the time lag data. We construct a corridor 
between two curves of time lag vs energy which includes all 
of the time lag data points (see Fig. |2). This allows us to 
constrain nf and nf. The lower and upper curves in Figure 
|2] correspond in panel (a) to nf = (6.1 - 8.0) x 10 18 cm 3 
and nf = (1.6 - 2.6) x 10 18 cm" 3 , respectively. For the 
sources in panels (b) XTE J1751-305 and (c) SAX J1808.4- 
3658 the density variations are nf = (6.0 - 6.6) x 10 19 cm" 3 
and nf = (2.9 - 3.6) x 10 19 cm" 3 . Thus, the plasma density 
of the "cold" reflector can change as much as 10 % during the 
entire data collection. 

4. PULSED FRACTION OF X-RAYS: OBSERVATION AND THEORY 

In Figure [3] we show the energy dependent pulsed frac- 
tion of (a) IGR J0029 1+5934, (b) XTE J 175 1-305 and (c) 
SAX Jl 808.4-3658. The pulsed fraction of IGR J0029 1 +5934 
vs energy is explained by the energy d ependent electron c ross- 
section cr e (E) = cr T (l - 2z) (see e.g. lPomranningiri973h and 
consequently by Compton cloud optical depth as a function of 
energy 

t c1 (£) = tt,ci(1-2z). (7) 

We assume that the accretion column is embedded in 
the Compton cloud as shown in Figure Q] (see also 
iTitarchuk et al. 1 120021 for the geometry details). Because 
T c i{E) decreases with energy, a larger fraction of the pulsed 
direct hard X-rays which originated in the accretion column 
can escape to the observer. This fraction is given by 

^rms,es(^') — 

A rms , es (0) e\p[-T c i(E)] (8) 

A different scenario of energy dependent amplitude forma- 
tion can be if there is no electron (Compton) cloud between 
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the accretion column (where the Comptonization spectrum 
is formed) and the observer. In this case the energy depen- 
dence of the amplitude is formed as a result of upscattering. 
The soft photons upscattered off hot electrons (in the accre- 
tion column) increase their energy with the number of scat- 
terings. On the other hand, the amplitude of the pulsed radi- 
ation exponentially decreases with the number of scatterings 
^nns(M) = A rms (0) exp{-/3u), where /3 is the inverse of the av- 
erage number of scatterings (see e.g. Sun vaev & Titarchukl 
1980) and consequently with an energy of the upscattered 
photon (see Eq. 



1 rms,up 



[«(£)] = A rms , up (E s )(E/E s y 



(9) 



(where a c \ = /?/ '40 is approximately the energy index of the 
Comptonization spectrum). A characteristic seed (disk and 
NS) photon energy is E s < 3 keV. Consequently, all photons 
at higher energies (> 3 keV) are produced by upscattering, 
and their amplitude should decrease with energy. 

In the general case, some of the pulsed soft photons, (a frac- 
tion A lms up ), may upscatter off hot electrons in the accretion 
column on the way out as other pulsed photons (a fraction 
Arms.es) forming the hard X-ray tail may escape to the observer 
while passing through the Compton cloud. In this case, the 
following combination of Eqs. ([8]and|9]l 



ArmsC^) — ^rms,up( E j E ^) 



+ A g-^i® 

^ ^H"ms,es c 



(10) 



lead to the formula of the emergent pulsed amplitude. 

In Figure|3]we present the best-fit model (see Eq. ITOb along 
with the data. We fit only those data points which correspond 
to energies higher than the seed photon energy. We find for 
(a) IGR J0029 1+5934 r T , c i = 3+0.2 and A rms>es = 98^ %, 
Arms.up = 5.2 + 0.9 %. For (b) r T , c i = 0.25 + 0.09 and A ms , £S = 
3.5 + 0.8 %, A^ up = 1.8 ± 1.1 % and (c) t t c] = 0.4 + 0.2 and 



5.9 + 6 



Mms.up 

all of the fits a c \ was ~ 0.8. 



2.5 ±1.2%, respectively. For 



5. DISCUSSION AND CONCLUSIONS 

In this Letter we have shown that the short time lags, phase 
shifts of the order of 100 /is between X-ray pulses at different 
energy ranges, detected in the accreting MSPs are the result 
of photon upscattering and downscattering in the plasma en- 
vironment of the NS. In our model the Comptonization emis- 
sion area, NS surface, and accretion disk are responsible for 
the time lag formation. The photon time lag correlations as 
a function of the electron densities of the reflector, n™ f , and 
the accretion column, n^ ', at a given energy are successfully 

reproduced for a reasonable choice of model (Cm, kT" ) and 
best-fit photon spectral parameters (a, kT 1 ™'). As illustrated 
in the schematic view of the accretion flow (see Fig. [TJ, the 
hard X-ray photons are downscattered in the relatively cold 



plasma of the disk and therefore produce the observed soft 
lags. The observed hard time lags are a result of thermal 
upscattering of soft photons coming from the disk and NS 
off hot electrons of the Comptonized region. This conclusion 
is based on the spectral shape of the emergent spectrum and 
time lag behavior as a function of photon energy for three 
MSP sources. Indeed, we find for IGR J0029 1+5934 with a 
btmax ~ 100 jus that rff ~3x n\°\ indicating that the densi- 
ties of the relatively "cold" disk plasma and of the accretion 
column are on the same order of 10 18 cm~ 3 for Cm ~ 0.1. For 
SAX J1808.4-3658 and XTE J1751-305 the observed energy- 
dependent time lags need disk density of a factor of seven 
higher. The hard lags of the emission due to Compton upscat- 
tering are negligible because of the expected compactness of 
the accretion column region. 

However, one can notice in Figure|2]that there are signs of 
turnover (saturation) for time lags vs energy at ~20 keV for 
SAX J1808.4-3658 and XTE J1751-305 which are similar to 
that in IGR J0029 1+5934 which has a turnover at about 7 ke V. 
The high values of n^ ' with respect to nf { would not allow us 
to see the contribution of the hard lags in the resulting time 
lag vs energy plot. 

It is also important to emphasize that the energy-dependent 
time lags allow us to reveal the plasma density variation 
within 10% near the NS. It implies that the mass accretion 
rate and ultimately the observed luminosity should also vary 
within 10%. Thus, we can conclude that the observed vari- 
ability of X-ray radiation in the compact objects is presum- 
ably driven by this variation of the density of the mass accre- 
tion there. 

Measureme nts of time lags have be e n done for black hole 
binaries (e.g. iKvlafis & Klimisl Il987t iMTUerl [l99l iNowakl 
fl995l) . Making a comparison between these measurements 
and that of time/phase shifts between X-ray pulsations in ac- 
creting MSPs is difficult. In the case of black hole binaries, 
the time lag is generally calculated in two different energy 
bands and is related to different Fourier frequencies. For black 
holes, the Fourier frequency dependence can be due to hydro- 
dynamical disk instability, mass accretion configuration, or 
quasi-periodic instability. This work with accreting millisec- 
ond NSs provides a unique advantage since the X-ray pulsa- 
tions give a perfect clock which allows to determine precise 
information about the photon delay at different energy bands. 
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